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Abstract 



We use a semiclassical approximation to derive the partition function for 
an arbitrary potential in one-dimensional Quantum Statistical Mechanics, 
which we view as an example of finite temperature scalar Field Theory at 
a point. We rely on Catastrophe Theory to analyze the pattern of extrema of 
the corresponding path-integral. We exhibit the propagator in the background 
of the different extrema and use it to compute the fluctuation determinant 
and to develop a (nonperturbative) semiclassical expansion which allows for 
the calculation of correlation functions. We discuss the examples of the single 
and double-well quartic anharmonic oscillators, and the implications of our 
results for higher dimensions. 
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I. INTRODUCTION 



The development of nonperturbative methods for quantum field theories remains as much 
a challenge as a necessity. The many areas of Physics where field-theoretic approaches have 
proved useful are full of examples of phenomena whose description lies outside the scope of 
perturbation theory. In most reliable treatment of the nonperturbative dynamics is 

still lacking. 

With the exception of a limited number of exactly solvable models, most systems re- 
quiring nonperturbative methods have been the object of approximate treatments. This is 
an unavoidable consequence of the complexity of the dynamics involved. Nevertheless, it 
should not prevent us from having a clearer understanding of the physics, as long as we have 
some control over the approximations made. 

The approximate methods most commonly used could, perhaps, be classified into three 
major groups: expansions around special limits, variational methods and numerical tech- 
niques. Semiclassical and expansions, Gaussian and Hartree-Fock methods, and Monte 
Carlo techniques in lattice calculations exemplify each respective group. In all of them, 
serious technical difficulties stand in the way of establishing reliable controls and, in some 
cases, of even fully implementing the approximation scheme. 

Semiclassical expansions certainly pose hard problems for both their implementation and 
control. In this article, we discuss a physical situation where such problems can be overcome: 
one-dimensional Quantum Statistical Mechanics or, from a field-theoretic point of view, 
scalar finite temperature Field Theory in zero spatial dimension. The drastic reduction in 
dimensionality is the price we have to pay to accomplish our goal: a thorough semiclassical 
treatment whose results can be compared to those of perturbation theory, as well as to 
exact ones, when available. One should not forget, however, that numerous applications in 
Physics can be studied within this framework, and that, as we shall indicate in the sequel, 
symmetries may reduce higher dimensional problems to situations where our results may 
still be used. 

The article will start with a discussion of the quadratic semiclassical approximation (Sec- 
tion IlQ), comprising an analysis of the extrema of the classical action and of the fiuctuation 
determinants around them. A connection to Catastrophe Theory allows for the identification 
of the minima amidst the complexity of the set of classical solutions. The next step (Section 

III) is the derivation, using techniques described in of the semiclassical propagator in 
the background of any one of the extrema, a crucial result which leads to the expression for 
the fiuctuation determinant and serves as the basis for a semiclassical expansion (Section 

IV) . The expansion provides, already in first order, nonperturbative estimates for correlation 
functions. Ground states of quartic anharmonic oscillators can be computed to that order: 
the result contains a sum over an infinite number of diagrams in perturbation theory, an 
indication of the power of the technique. Finally (Section V), we summarize our conclusions 
and comment on possible applications of these results in higher dimensions. 



This section is strongly based on Ref. [||]. 
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II. THE QUADRATIC SEMICLASSICAL APPROXIMATION 



A. Extrema and determinants 

The partition function for a one-dimensional quantum-mechanical system in an arbitrary 
potential, V{x), can be written path-integral 

foo i-x{l3h)=xo 

Z{(3)= dxo [DxiT)]e-'/\ (2.1) 

J —oo J x(0)=xo 



r/3h 

S= dr 
Jo 



^Mx^ + V{x] 



(2.2) 



As mentioned before, this can be viewed as the partition function for a scalar finite temper- 
ature field theory, in the limiting case where the number of space dimensions is zero. 

We may expand the Euclidean action around its minima, i.e., the stable classical trajec- 
tories,0 and keep terms only up to quadratic. The functional integral over the paths can, 
then, be performed to yield a quadratic semiclassical evaluation of (^7^ 

N{xo,f3) 

E 



/•oo ^(^0,/i<) 

Zsc{P)= dxo E e-^^(^-^)/^A7^/'(a;o,/3), (2.3) 

J —oo ■ 1 



where Sj{xo,P) denotes the action and Aj{xo,P) represents the determinant of the fluctua- 
tion operator, 

F, = -M^^+V"[x,{r)], (2.4) 

both calculated at the j-th classical trajectory, xji^r), satisfying the boundary conditions 
Xj{0) = Xj{Ph) = Xo] N{xo,P) is the number of classical trajectories which are minima of 
the action functional. 

The action has a simple expression in terms of the turning points of the classical trajec- 
tory: 



± 

XQ 



Sj{xo,f3) = (3hV{xi)±2 dxMv{x,xi) + 2n . dxMv{x,xi). (2.5) 



Here, v{x,y) = ^J{2/M)[V{x) — V{y)] and x-^ (xl) are turning points to the right (left) of 
xq. The first term in (|2.5|) corresponds to the high-temperature limit of Z{f3), where the 
classical paths collapse to a point, i.e., Xq. For motion in regions where the inverted 

potential is unbounded (hereafter called unbounded motion), n = 0. For periodic motion, n 
counts the number of periods and the second term is absent. For bounded aperiodic motion, 
all terms are present. The last two terms will be negligible for potentials which vary little 



^For classical trajectories we mean solutions of the Euler-Lagrange equation Mx — V'{x) = 0, 
which is the equation of motion for a particle moving in the potential minus V. 
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over a thermal wavelength, A = j3/M. However, by decreasing the temperature, quantum 
effects become important. 

For trajectories having a single turning point (n = 0), are given implicitly by 



Jx 

and the ffuctuation determinant by 



dx 



(2.6) 



Aj{xo,P) 



Anh^lVixi) - V{xo)] 



MV'ix^ 



8(3 
dxi. 



(2.7) 



Xo 



For trajectories with two turning points, we may generalize (|2.7| ) to 



A,{xo,P) 



A7rh^[V{xi) -V{xo)] J 



M 



V'{xi 



d/3 
dx\ 



n_ 



-I xo 



V'ixi 



dl3 
dx''_ 



(2.8) 



n± counts the number of visits to the two turning points, a;±. As we shall see in the sequel, 
trajectories with more than one turning point are naturally excluded from the calculation. 
The derivations of ( ^.7] ) and ( |2.8| ) will be outlined in section III. 

Formulae ( |2.5| ), (|2.7| ) and (|2.8|) do not require full knowledge of the classical trajectory, as 
the dependence on Xj{T) comes only through the turning point. Furthermore, equations ( p.7|) 
and (|2.8| ) yield a direct evaluation of the determinant which bypasses solving an eigenvalue 
problem for each stable classical trajectory. To see how this can simplify the evaluation 
of Zsc(/3), let us take the harmonic oscillator, V{x) = \Muj'^ example. In this 

case, given xq and (3 there is only one trajectory, with a single turning point, given by 
x+(x_) = Xo/ cosh(/3^/2) for xq < (> 0). 5" and A can also be readily calculated, the 
final result being 



dXQ e'^^'^^o/^^ tanh(/3?iLj/2) 



(2.9) 



which in this case is exact. 

For the single- well quartic anharmonic oscillator, V{x) — ^^.^^ ^ , ^ 
still true that there is only one trajectory, given xq and /3, which has a single turning point. 
Eq. (P;BI) yields: 



, A > 0, it is 



x± = Xo cn(M-|-, /c-t 



'1 + 



Ax^ 
Mcj2 



k. 



2 _ Mu'^ + (l/2)Ax| 
+ Ax^ ' 



(2.10) 
(2.11) 

(2.12) 



where cn(M, /c) is one of the Jacobian elliptic functions 0. We may, then, use Eqs. ( p.5| ) 
and ( |2.7| ), and change the integration variable in ( pl3|) from xq to x± to obtain Z^dP)- The 
resulting integral can be evaluated numerically. 
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A more interesting situation occurs in the case of the double-well quartic anharmonic 
oscillator, V{x) = A(x^ — a^)^. For > a^, only single paths with single turning points 
exist for fixed Xq and f3. However, there is also a region, < a^, where the classical motion 
is bounded and a much richer structure exists , one in which more than one classical path 
may exist for given values of xq and (5. 

In a region of bounded classical motion (a well in — V), such as < for the anharmonic 
oscillator, the number of classical trajectories changes as the temperature drops. If < < 
n/hum (where cj^ = \J —V"{xm) / M and Xm is a local minimum of — V^), for every Xq in 
this region there is only one closed path, with a single turning point, satisfying the classical 
equations of motion. For xq < Xm (> a^m) this path goes to the left (right) and returns to 
Xq. For Xq = Xm, it sits Still at the bottom of the well. It is this single-path regime which 
goes smoothly into the high-temperature limit. 

For P = Tr/hujmi the solution that sits still at Xm becomes unstable. Its fluctuation 
operator is that of a harmonic oscillator with uj'^ = uj"^. Its finite temperature determinant 
is A{xm, P) = 2'7ihsm{Phum)/Mum- This goes through zero at (3 = n/hum and becomes 
negative for P > Tr/hum, thus signaling that x(r) = Xm becomes unstable. At the same 
time, two new classical paths appear, symmetric with respect to Xm, as depicted in Fig. 1. 
Therefore, at Xq = x^, we go from a single-path regime to a triple-path regime as we cross 
P = n/hum- The two new paths are degenerate minima, whereas the path that sits still at 
Xm becomes a saddle-point of the action, with a single negative mode. 

As P grows beyond ir/hum, an analogous situation occurs for other values of xq inside the 
well. When the fluctuation determinant around the classical path for a given xq vanishes, 
a new classical path appears. Its single turning point lies opposite, with respect to Xm, 
to that of the formerly unique path. If P is increased further, this trajectory splits into 
two, as illustrated in Fig. 2. One is a local minimum of the action, whereas the other is a 
saddle-point, with only one negative mode. Again, we have transitioned from a single to a 
triple-path regime. As P grows, the triple-path region spreads out around Xm- The frontiers 
of that region are defined by the points xq such that {dp / dx±)^^^ = 0, where the expression 
( p.7| ) for the determinant vanishes. 

The phenomenon just described is an example of catastrophe. It takes place whenever the 
number of classical trajectories changes, with one or more of them becoming unstable, as we 
lower the temperature. Conversely, we may say that the phenomenon is characterized by the 
coalescence of two or more of the classical trajectories, as we increase the temperature. This 
is akin to the ocurrence of caustics in Optics where light rays play the role of classical 
trajectories and the action is replaced with the optical distance. In the next subsection, we 
shall use the mathematics of catastrophes to understand how the set of classical trajectories 
changes with temperature. 

B. Catastrophe Theory 

If we use the basis of eigenfunctions of the fluctuation operator around classical paths, 
and denote by si the coordinate associated with the direction of instability of our example, 
the action can be viewed as S* = S{si, . . . ; xq, P), with the dots referring to all others. 
Catastrophe Theory allows us to write down the "normal form", S^, of this generating 
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function in the three-dimensional subspace made up by the unstable direction of state space 
(the set of paths) and the two control variables, f3 and Xq] it is this subspace which is relevant 
for the study of the onset of instability. As catastrophes are classified by their codimension 
in control space, which here is two-dimensional, only those of codimension 1 (the fold) or 2 
(the cusp) are generic (i.e., structurally stable). The pattern of extrema then leads to the 
cusp, whose generating function is 

Sn{si; Ml, vi) = sl + uisl + visi , (2.13) 

where ui and vi, the control parameters, are related to (3 and xq, respectively. This catas- 
trophe is defined by 



dSN d^SN 



dsi dsj 



0. (2.14) 



Eliminating si from these equations, we can draw the bifurcation set in control space (the 
curve vf = — 4uf/27), as well as the pattern of extrema of Sn (Fig- 3). Assuming that ui and 
Vi are linearly related to {(3 — n/hum) and (xq — Xm), respectively, for small values of both, 
we can also schematically plot the classical action (i.e., the action for classical trajectories) 
as a function of Xq for different values of see Fig. 4(a,b). 

New classical trajectories, accompanied by new catastrophes, will occur as we keep in- 
creasing p. From f3 = 2'K/huJm, we start having trajectories with x{0) = 0, i.e., with one full 
period. For these, the determinant of fluctuations vanishes because the first, not the second 
bracket of ( ^.Tj ) goes trough zero. However, the fluctuation operator around such trajecto- 
ries already has a negative eigenmode. This follows from the fact that the zero-mode, given 
by Xci{t), has a node. Thus, it cannot be the ground-state for the associated Schrodinger 
problem. This means that a new catastrophe takes place along a second direction in function 
space. 

This new catastrophe is also a cusp. At Xm, as P becomes larger than 2n /hum, the 
solution that sits still acquires a second negative eigenvalue — thus becoming unstable 
along a new direction in function space — and two other solutions appear. They both 
have a negative eigenvalue along the first direction of instability and a positive eigenvalue 
along the new direction of instability. If we combine the information from both directions, 
we find that we go from two minima and a saddle-point with only one negative eigenvalue 
(hereafter, a one-saddle) to two minima, two one-saddles and one two-saddle. The two one- 
saddles are degenerate in action and represent time-reversed periodic paths. As f3 increases 



^This plot can be obtained by exploiting the relation between the cusp and the swallow- 
tail catastrophes. The latter has a generating function whose normal form is W{si;a,b,c) = 
sf + asi + bsi + csi, where a, b, c are control variables. The extremum condition is, then, 
dW/dsi = 0. The identifications Sjy = — c/5, ui = 3a/5 and vi = 2b/5 will make the additional 
condition that defines the swallowtail, d'^W/dsl = 0, coincide with the requirement dS]\f/dsi = 0. 
(In the jargon of catastrophe theory: the bifurcation set for the swallowtail coincides with the 
equilibrium hyper surf ace of the cusp.) 
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beyond I-k jhuora-, the same phenomenon takes place for xq around Xm- the one-saddles that 
already existed for /? < 2ti /hujm in the three-path region around Xm become unstable along 
a second direction, and a five-path region grows around x^, with two minima, two (periodic, 
degenerate in action) one-saddles and a two-saddle. 

The new cusp can be cast into normal form using the second direction of instability: 

Sn{s2; M2, V2 = Q)=sl + U^sl ~ {Sl + y )2 . (2.15) 

The absence of the linear term in Eq. ( p.l5| ) [compare with Eq. ( p.l3| )] reflects the degeneracy 
in action of the two one-saddles. However, there is another degeneracy: since V does not 
depend on /5, if Xc\{t) is a solution of the Euler-Lagrange equation, so is Xc\{t + Tq). If 
Xc\{t) is periodic, Xc\{t + Tq) describes the same path (and so has the same action), but 
with another starting point. This can be represented in Eq. ( [^.15D by the choice U2 oc 
[(xo - Xmf + k{2'n/hujm - 13)], with A; > 0. See Fig. 4(c) 

As we approach (3 = Sir/hujm, the situation becomes similar to the one near f3 = Tr/huOrn- 
The difference is that we now have to deal with a path with more than one and less than 
two periods. Near f3 = An/tujm, two-period paths intervene, a situation similar to that at 
P = 2t[ /hwm- The pattern which develops [1^ is depicted in Fig. 5. 

Despite the fact that we keep adding new extrema as we lower the temperature, only 
two of these are stable (i.e., minima). In a semiclassical approximation in euclidean time, 
appropriate to equilibrium situations, these are the only extrema we have to sum over, 
meaning that N{xo,P) never exceeds two. For /3 > n/hum, there will be regions with either 
one or two minima of the classical action. The transition from a single-minimum to a double- 
minimum region occurs at values of xq where the fluctuation determinant vanishes due to 
the appearance of a caustic, thus leading to a singularity in the integrand of Eq. (p.3|). This 
is not a disaster, however, as this singularity is integrable. (Such a singularity is an artifact 
of the semiclassical approximation. It disappears in a more refined approximation PJTT|-p!3|, 
in which one includes higher order fluctuations in the direction(s) of function space where 
the instability sets in.) We shall exhibit the details of this calculation in a forthcoming 
pubhcation p^ . 

To conclude, two remarks: (i) our results may be used to derive a semiclassically improved 
perturbation expansion, which will be presented in Section IV; (ii) the semiclassical partition 
function incorporates, in an approximate way, quantum effects that become more and more 
relevant with the increase of the thermal wavelength. At very low temperatures, however, 
the quadratic approximations inherent to the semiclassical approximation fail to capture 
the details of the potential, which are important in the regime of large thermal wavelengths. 
Thus, we expect our calculations to be a better approximation at high temperatures. Indeed, 
studying the T = limit of Z^dli) for the anharmonic oscillator, one finds no corrections 
to the unperturbed ground state energy. Nevertheless, this can be corrected by using the 
semiclassically improved perturbation theory just mentioned. 



III. THE SEMICLASSICAL PROPAGATOR 

In the next section we will show that, in order to go beyond the quadratic approximation 
of the preceding section, it is necessary to derive the semiclassical propagator, defined as: 
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Gj{vuV2;n,T2) = r^"'^ ''[Dr/(r)]e-^^/'* , (3.1) 
"''7(n)='?i 

2 

,2 



(3.2) 



Sj involves the fluctuation operator, Fj, which corresponds to the second functional deriva- 
tive of the action around the classical trajectory, Xj(r). 

Formula (|3.1|) is the path-integral expression for the Euclidean time evolution operator, 
Pj, of the quantum-mechanical problem defined by 

-h-^i,,(r,0) = Hj{T}p,{T.O), (3.3) 

Indeed, the propagator is given by 

Gj{vi,V2;ri,T2) = {v2\pj{T2,Ti)\rii) . (3.5) 

Thus, it is just a density matrix element, as suggested by equation (|3.4|) . 

As the integral in ( |3.1[ ) is Gaussian {Hj is quadratic), it is completely determined in 
terms of the extremum, f]j{T), of Zlj[?7(r)], which satisfies 

-M^ + V"[x,{r)]rj,=0, (3.6) 
with the boundary conditions 

Vjin) = Vi and f]j{r2) = r]2 . (3.7) 
The semiclassical propagator is, then, 

G,(r/i, ri2; n, ra) = G,(0, 0; n, T2) e"^^/^ , (3.8) 
with = T,j[fjj{T)]. Using equations ( p.l|) , (|3.6|) and ( ^^ , we have 



(3.9) 



The fluctuation determinant discussed in the previous section can be obtained as 

Aj'/\xo,P) = G,{0,0;0,(3h) . (3.10) 

The extremum, f]j{T), will be the linear combination of two linearly independent solutions 
of equation (^.6|), aji^r) and bj{T), which satisfies the boundary conditions ( |3.7| ): 

^,^r) = ^!^B^^aAr) + ^^^l^^^b,ir), (3.11) 
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with Qji = aj{Ti) and bji = bj{Ti). This allows us to compute Hj. Introducing the functions 



JF,i(r) = ajibj^r) - aj{T)hji , 
JF,-2(r) = aj2bjir) - aj{T)hj2 , 

which satisfy jFji(ri) = 0, J^j2(t2) = and J-'ji{T2) = — ^j2(Ti), and the Wronskian 

Wjir) = %(r)6,(r) - a,(r)6,(r) , 



we may write 



M 



a,b=l 



where 



„{i) — 3 



-^72(7-1 



-^72(^1; 



1 ^12 — 



Si) 



^i2(n) ' -^ii(^2; 

Using equations ( p.4| ), (p.3|), ( |3.5| ) and (|378|), we obtain 

_c'r2 2 ^ji(,r2jj 
which leads to 

G,(0,0;ri,r2) = Q(ri)[-F,i(r2r^/^ 
As r2 ^ Ti, we must recover the free-particle result: 



) ■'22 



•^il(^2) 



lim G,(0,0;ri,r2) 



M 



27rh{T2 - Ti) 



1/2 



(3.12) 
(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



Expanding ^^1(^2) in equation ( p.l8| ) around ti, and using equation ( p.l9| ), we can determine 
Cjln). Finally: 



G',(0,0;ri,r2) 



27rh J^jl{T2) 



1/2 



(3.20) 



All that remains is to find the two linearly independent solutions of (3^), cij{T) and bj{T). 
The arguments of section II allow us to restrict our attention to single turning point paths, 
although we will sketch how one should proceed in the most general case, at the end of this 
section. 

It can be easily checked that xj^t) satisfies equation (|3.6|), by simply differentiating the 
Euclidean classical equation of motion with respect to Euclidean time. Then: 



(3.21) 
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If we consider r < /3^/2, one can check that 

also satisfies equation ( p.6|) , again thanks to the equation of motion satisfied by Xj{T). 
Therefore, it can be used as bj^r) for r < i3h/2. It can be shown that c± = Cj{(3h/2) and 
c4 = Cj{Ph/2) are well-defined [l^. Clearly, these will be the values of bj and bj at [3h/2. 
In order to find our other solution, for r > /3h/2, we cannot use Cj(r), which is ill-defined 
beyond [3Ti/2. However, we may use 

A linear combination of aj(r) and Cj(r) may, then, be constructed whose value at (3h/2 and 
that of its derivative coincide with bj{Ph/2) = Cj{[3Ti/2) and bj{j3h/2) = Cj{ph/2). These 
two conditions determine the linear combination and guarantee continuity. As a result, we 
may now define &j(r) for all values of r: 

^^■^^)-|x,(r)/;,^-2c4c4x,(r) for r > f3n/2. ^^'^^^ 

We note that both aj (r) and bj (r) can be completely expressed in terms of the potential and 
of the value of the (single) turning point. Also, for the general case of two turning points, 
we can use the above procedure to obtain a piecewise construction of bj^r). For the first 
piece, we trade f3h/2 for the r* where the first turning point is visited. Then, (3h will be 
traded for 2t*. This will allow us to define &j(r) up to 2r*. We may proceed in the same 
fashion to the next turning point and repeat the operation to account for the number of 
cycles involved. 

Equation ( |2.8|) was obtained from Gj{0, 0; 0, l3h), for the general case. From the knowl- 
edge of f]j{T), we can write a final expression for Gj{rii,ri2; Ti,T2): 

Gj{vi,V2;ri,T2) = ^==e ^ (3.25) 
^2na] 

^. = ^•^.1(^2), (3.26) 

A{2 = ^^^^^^ + ^^^^^^^i^^^^^ (3 27) 

aj(r2) 

This will be the basis for the semiclassical expansion to be derived in our next section. 

IV. THE SEMICLASSICAL EXPANSION 

We are now in a position to go beyond the quadratic approximation. To do so, we expand 
the action around a classical trajectory, Xj^t), to all orders. Letting x(r) = Xj(r) -|- rj^r), 
with 77(0) = ri{Ph) = 0, we obtain: 
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S[xir)] = S[x,iT)] + - drr^{r)F,r^{T) + Y.- y(")[x,(r)] r/», (4.1) 

'J ^ 2 vt. JO 



where ^'•"^ = d"'V/dx'^, and the second term is the quadratic piece. Use was made of the 
fact that higher functional derivatives involve products of delta functions. 

As we take e~^^^ and integrate over the deviations, //(t), from the classical trajectory, 
we may expand the higher than quadratic terms in a power series. Defining: 



/ dT5L,^Y.-J rfrV(")[x,(r)]r^"(r), (4.2) 
JO n. JO 



we may interchange the integral over r above with the sums over deviations from the classical 
path. This has the effect of breaking up the path integral into pieces that go from zero to 
the various times corresponding to multiple insertions of 6Lj and, finally, to (3h. The result 
of doing the path integral in each piece yields the semiclassical propagator at different time 
intervals. The expression for the partition function becomes: 

Z{(3)= dxo E e-'^/^[G,{0,0;0,(3h) 

J-oo 

i'f3h poo 

dTi / dr]iGj{0,r]i;0,Ti)6Lj{Ti,rii)Gj{r]i,0;Ti,(3h) 



oo 



h Jo 

I i-lSh POO 



2h^ 



rpfi foo rT2 poo 

/ dT2 / dr]2 / dri / dr]iGj{0,rii;0,Ti) 5Lj{Ti,T]i) 

JO J —oo JO J —oo 



xG,{r],,r]2; n, r2) 6L,{t2, m) G,{r]2, 0; ts, ph)] + 0(3) (4.3) 

Clearly, the Gaussian dependence of the semiclassical propagator on the r] variables, together 
with the fact that SLj only involves powers of rj, allows the integrals over t] to be easily 
performed. As for the r dependence, it is completely known, although the integrals involved 
may not be done analytically. The integrals over r can, alternatively, be transformed into 
integrals over x (the classical position). We should note that the first term in ( |0| ) is just the 
fluctuation determinant and that each individual term in the series corresponds to an infinite 
sum of diagrams in ordinary perturbation theory, plus additional terms which depend on 
the classical solution. 

In order to illustrate the usefulness of the expansion, let us outline some of its conse- 
quences in the case of the single- well quartic anharmonic oscillator, V{x) = ^mcu'^x'^ + ^x^, 
already discussed in section II. If we define a dimensionless coupling g = \h/2m'^u^, the 



first order correction in the expansion, i.e., the second term in the bracket of equation ( [4. 3D , 
may be written as —3g Gj{0, 0; 0, (3h) Ij{l3, g), where Ij is a dimensionless integral: 

I, = j^^^, I dr [b%r) + 4c4c465(r)a,(r) + (2c4ci)^6|(r)a|(r)] (4.4) 

The functions aj{T) and bj{T) can be expressed in terms of Jacobian elliptic functions and 
elliptic integrals |]14| which we shall not write down. All we want to emphasize is that g) 
can be expanded in a power series in g, just as Gj{0,0]0, ph). In ordinary perturbation 
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theory, the first order term is of order g, and goes at most as f3, for large f3. An n*'* order 
term will go like (7"/?" (subleading corrections in f3 are present). Here, in first order, we have 
a whole series in g times (3. Thus, if we remember that the exact ground state EQ{g) is the 
only state to contribute at large /3: 



e 2 



Eo{g) - 



tiUJ 



+ 0{(3') 



(4.5) 



we see that we obtain a nonperturbative estimate for EQ{g), while perturbation theory would 
only give us an order g correction. 

In order to compute correlation functions, one might introduce an external current, j(t), 
coupled to the variable, x{t), and use the knowledge of the semiclassical propagator to 
integrate the linear plus quadratic piece of the action to derive a functional of j (r), to be 
used in the expansion in the standard way . 



V. CONCLUSIONS AND OUTLOOK 

The methods developed in the preceding sections allow us to deal with problems in one- 
dimensional Quantum Statistical Mechanics by means of a semiclassical expansion which is 
fully calculable and yields nonperturbative results. 

Problems in higher-dimensional Quantum Statistical Mechanics with potentials that have 
spherical symmetry may be reduced, by means of partial-wave decompositions, to one- 
dimensional problems wherein our methods could be applied. 

Furthermore, in finite temperature Field Theory, where the point xq of our previous treat- 
ment is to be replaced with the (static) field configuration, (j){x), we already see the need to 
consider extrema which must depend on r, the Euclidean time, as we lower the temperature. 
Catastrophes will certainly occur, forcing us to rethink commonly used approximations that 
use constant field configurations as backgrounds to derive effective potentials. Here, again, 
under very special circumstances, one might hope that the methods developed so far may 
shed some light into the intrincacies of nonperturbative phenomena. 
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Figure Captions: 

Figure 1: (a) Classical paths at Xm for /3 < ■n/fiujm (1) and /3 > T:/fiuJm (2 and 3); (b) 
Sketch of how the extrema change along the unstable direction in function space. 

Figure 2: (a) Classical paths at x„i for (3 < (3c (1), (3 = (3c {1 and 2&3) and (3 > (3c (1,2 
and 3), n/hum < (3c < 2-k /fujjm'i (b) Sketch of how the extrema change along the unstable 
direction in function space. 

Figure 3: Bifurcation set for the cusp; pattern of extrema shown schematically. 

Figure 4: Evolution of the classical action as (3 changes, (a) < /9 < 'n/TyjJm'i (b) T:/fiu)m < 
(3 < 2t{ /fujJra') (c) 2'n/huim < (3 < STr/huJm- 

Figure 5: Partition of control space into p-solution regions (p = 1, 3, 5, 7, . . .). 
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